Resonance fluorescence in driven quantum dots: electron and photon correlations 
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We study the counting statistics for electrons and photons being emitted from a driven two 
level quantum dot. Our technique allows us to calculate their mutual correlations as well. We 
study different transport configurations by tuning the chemical potential of one of the leads to find 
that the electronic and photonic fluctuations can be externally manipulated by tuning the AC and 
transport parameters. We also propose special configurations where electron-photon correlation is 
maximal meaning that spontaneous emission of photons with a well defined energy is regulated by 
single electron tunneling. Interesting features are also obtained for energy dependent tunneling. 
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I. INTRODUCTION 

Quantum dots offer an ideal playground for testing co- 
herent and quantum optical effects in an artificially de- 
signed solid-state environment^, with the additional ben- 
efit of having electronic transporlj^ii^ as a 'spectroscopy' 
by coupling to external fermionic reservoirs and counting 
the flow of single electrons. Complex behaviour emerges 
through electron-electron interactions and the interaction 
between electrons and other bosonic excitations such as 
photons* or phonons^. 

The exchange of ideas between mesoscopic transport 
and quantum optics has turned out to be quite fruit- 
ful. For instance, thermal electron antibunching was 
observed experimentally by performing flanbury Brown- 
Twiss-type experiments in mesoscopic conductors^. This 
fermionic antibunching has also been used for regular 
photon sources in p-n junctions^ and quantum dots^'». 
Reversely, bosonic statistics can be studied in quantum 
conductors such as beam splittersiS, nanoelectromechan- 
ical systemsii or quantum point contactaiS,, where pho- 
ton antibunching was predicted^^. Another example is 
coherent population trapping and dark states in multi- 
level atoms^ii^, originally proposed for driven three level 
quantum dotsi^ and then extended to triple quantum 
dots in a simple triangle configurationii. 

Quantum transport also benefits from the adoption of 
theoretical tools that are well established in quantum op- 
tics. Specially relevant in the last years has been the de- 
velopment of noisei^ii^ and full counting statistic a^°'^^i^^ 
for electrons. Here, many of the relevant ideas and 
techniques were in fact originally developed in quantum 
optics in the context of counting single photons that 
are emitted from a single atomic sourc o^'^'^'^i^^'^^i^^'^^ . 
Recently, electronic counting statistics has become ex- 
perimentally accessible for incoherent transport through 
quantum dot (QD) systems by the analysis of the time- 
resolved current flowing through a quantum point con- 
tact electrostatically coupled to the m^^'^^i'^^'^^ . ffowever, 
the backaction of the quantum point contact on the QD 
destroys its internal coherence. Though traces of coher- 



ence have been measured in shot noise through stacks of 
double quantum dots^, the access to higher order cum- 
mulants is still a challenging problem. 

Our aim in this work is to study the influence of elec- 
tronic transport on the photonic emission statistics in a 
quantum dot system, and vice-versa. Two-level systems 
give particularly interesting features both for optical and 
transport quantities: in optics, resonance fluorescence in 
two level atoms is the simplest case of a quantum photon 
source where photon antibunching occurs^, fn trans- 
port, quantum dots with two or more capacitively cou- 
pled levels show electronic bunching in dynamical chan- 
nel blockade configuration^. As will be shown, these 
properties can be studied in a two level quantum dot 
which is illuminated by a resonant ac electric field, where 
bosonic resonance fluorescence (due to phonon or pho- 
ton mediated relaxation processes) is modified by elec- 
tronic transport, and dynamical channel blockade de- 
pends on both coupling to a boson bath'^^ and the driving 
parameters'^^. 

We show how the combined statistics of Fermions and 
Bosons is a very sensitive tool for extracting informa- 
tion from time-dependent, driven systems. In particular, 
phonon emission has been measured by its influence on 
the electronic current in two level systems^. We analyze 
the electron and photon noises and find that they can be 
tuned back and forth between sub- and super-Poissonian 
character by using the strength of an ac driving field or 
the bias voltage. 

For this purpose, we develop a general method to si- 
multaneously extract the full counting statistics of sin- 
gle electron tunnelling and (boson mediated) relaxation 
events, as well as their mutual correlations. 

This paper is organized as follows: In Sec. |TT] we 
present our system and how to obtain the counting statis- 
tics of electrons and photons which are calculated in 
Sees, mil IIVI and |V] for different chemical potentials in 
the right lead. In Sec. IVI[ we present a special con- 
figuration where electron-boson correlation is maximal. 
The energy-dependent tunneling case is studied in sec- 
tion IVIII Finally, Sec. IVIIII presents our conclusions. 
Due to the length of some of the analytical results, we 



include some appendices with the coefficients that allow 
their calculation. 



II. MODEL AND TECHNIQUE 

Our system consists of a two level quantum dot (QD) 
connected to two fermionic leads by tunnel barriers. The 
Coulomb repulsion inside the QD is assumed to be so 
large that only single occupation is allowed (Coulomb 
blockade regime). The lattice vibrations induce, at low 
temperatures, inelastic transitions from the upper to the 
lower state. In analogy to resonance fluorescence in quan- 
tum optics, a time-dependent ac field with a frequency uj 
drives the transition between the two levels ei, £2 close 
to resonance, A^ = 62 — Si — uj ^ (we will consider 
e = h = 1 through the whole text), which allows us to 
assume the rotating wave approximation. Thus, the elec- 
tron in the QD is coherently delocalized between both 
levels performing photon- assisted Rabi oscillations^^. 
For simplicity, we consider spinless electrons. Studies 
of spin currents in a quantum dot coupled to a quantized 
driving or an AC magnetic field can be found in Refs. — 
and—, respectively. 

This system is modelled by the Hamiltonian: 
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where hqi, , Cka and di are annihilation operators of bosons 
(with momentum q and polarization v) and electrons in 
the lead a and in the level i of the QD, respectively, 
and O is the Rabi frequency, which is proportional to the 
intensity of the ac field. 

The terms proportional to Vai and A^-'^ in Eq. ^ give 
the coupling of the electrons in the QD to the fermionic 
leads and their interaction with the bosonic bath, respec- 
tively. In the following, we assume a basis of quantum 
dot levels where diagonal matrix elements of A^{, play a 
minor role and we set Xi^, = 0. For coupling to photons, 
this would correspond to a usual dipole interaction with 
the electromagnetic field. For phonons, this is justified if 
we are mainly interested in weak coupling and relaxation 
processes by spontaneously emitted bosons in the relax- 
ation from the upper to the lower level. The electron- 
boson coupling term, in the rotating wave approxima- 
tion, can then be written as X^qj/ ^qi' [dldiaqu + H.c. j. 
For the sake of illustrating our results with concrete phys- 
ical processes we will refer to photon emission in the fol- 
lowing, i.e. the bosonic bath corresponds to the photon 
vacuum without additional coupling to phonons. 

Finally, the coupling to the fermionic and bosonic 
baths terms are responsible for the incoherent dynam- 



ics and they can be considered apart in the derivation of 
the master equation for the reduced density matrix. 

Considering the quantum-jump approach^^ii^ to elec- 
tronic transport and boson emission events, the equation 
of motion of the reduced quantum dot density matrix can 
be written as 



p{t) = Co{t)p{t) + C,{t)p{t) + Cp(t)p{t), 



(2) 



where Ce and Cp are the Liouvillian jum,p superopera- 
tors responsible for the incoherent events: electron tun- 
neling from the system to the collector and relaxation 
by spontaneous photon emission. Thus, we can consider 
a density matrix resolved in the number of accumulated 
electrons in the collector, rig, and the number of emitted 
photons, ripi 



p(t)= 5] p("-">')(t). 



(3) 



where p("<!'"p)(t) gives the probability that, during a cer- 
tain time interval t, n^ electrons have tunneled out of a 
given electron-photon system and Up photons have been 
emitted. 

This allows us to define the generating function^^: 



G(t, Se, Sp) — 
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by introducing the electron (photon) counting variables, 
Se(p). The derivatives of G(t, Se, Sp) with respect to the 
counting variables give the correlations: 



dsids^ 
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(5) 
up to any order. 

One can derive the equations of motion for the gener- 
ating function as previously done for the density matrix: 



G(t, Se, Sp) — M (Se, SpjCrit, Sg, Sp), 



(6) 



that generalizes the master equation, p{t) — M(l, l)p{t). 
The jump superoperators affect only the diagonal ele- 
ments of the generating function in the same way that 
rate equations involve only the occupation probabilities 
-given by the diagonal elements of the density matrix. 
The electronic one acts as: 



{C,it)G{t,Se,Sp))^^ 
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where F^^ is the rate for processes that increase/decrease 
the number of electrons in the collector by a transition 
from state |fc) to state |m). The same can be done for 
photons, with the difference that the number of detected 
photons can only increase. Then, one only has to intro- 
duce the corresponding counting variables in those terms 
corresponding to the tunneling of an electron to the col- 
lector lead, and the emission of a photon. The relevant 



elements of the density matrix can be written as a vector, 
P = {Poo, Pii,Pi2, P2i,P22V, where poo gives the occupa- 
tion of the empty state, pn and p22 correspond to the 
ground and excited electronic states, respectively, and 
Pi2 and P21 — /3*2 are the coherences. Then, for the case 



where the tunneling barriers are equal for both energy 
levels, i.e. V^i = Va2, the equation of motion of the 
generating function ^ is described, in the Born-Markov 
approximation, by the matrijs^i^: 
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where, by further considering that the density of states in 
the leads is rather constant so da{u>mn) = da, the tunnel- 
ing rates through the lelft(right) lead are equal to rL(R) = 
27rdL(R)|^L(R)P- We will consider a high bias configura- 
tion where the chemical potential of the left lead is well 
above the energy levels and the transitions between the 
right lead (with a chemical potential /i) and the state i in 
the QD are weighted by the Fermi distribution functions 
f, = f{e,-fi) = (l-He(^'-'^)'5)"' and /, = !-/,. 7 is the 
spontaneous emission rate due to the coupling with the 
photon bath: 7 = 2ttY,,, J d^qg{q)\Xq„\'^S{\q\v ~ 82 + ei), 
where g{q) is the density of states—. The decoherence 
is given by A12 = -i ((A -I- /2)rR + 7). The Fermi en- 
ergy of the left lead is considered high enough that no 
electrons can tunnel from the QD to the left lead. All 
the parameters in these equations, except the sample- 
depending coupling to the photon bath, can be externally 
manipulated. 

Taking the Laplace transform of the generating func- 
tion, G{z, Se, Sp) — {z ~ M)^^p{0), where p(0) is the ini- 
tial state, the long-time behaviour is given by the residue 
for the pole near z — 0. From the Taylor expansion of 
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the pole zq = Yl 

write trG{t,Se,Sp) r^ g{se,Sp)e^°*' and obtain, from ([5]), 

the mean value, (n^i^p)), as well as the higher order cu- 



mulants, k 
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+ (CIO(OI) + 6C20(02) + 6c30(03))i, 



which give the variance and skewness of the probability 
distribution, respectively. In the large time asymptotic 
limit, all the information is included in the coefficients 



Cmn- Thus, we obtain the stationary current and the low 
frequency noise: 
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(10) 
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respectively. Then, the Fano factor is i^e(p) = 1 + 
2c2o(o2)/cio(oi) so that, the sign of the second term in 
the right hand side defines the sub- {F < 1) or super- 
{F > 1) Poissonian character of the noise. 



In the limit F 



L{R) 



0, the pure Resonance Fluores- 



cence case for the noise of the emitted photons, formally 
equivalent to the expression for emitted photons in quan- 
tum optics^® is obtained: 



Fp{r, = o) = i- 



2»2(37^-4Ag,) 
(72 + 2Q^ + 4A2 )2 



(12) 



yielding the famous sub-Poissonian noise result at reso- 
nance (A^j = 0). The detuning between the frequency of 
the AC field and the levels energy separation, A^, con- 
tributes to restore super-Poissonian statistics, as seen in 
Fig. [21 In the following, only the resonant case will be 
considered unless the opposite were indicated. 
Electron-photon correlations are obtained from: 

, %(1,1) %(1,1) 2 

{rLerip) = — — coii H — ci^t + cut + ciocoii 
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Then, by defining af^ = {uiUj) - {ni){nj), 
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(13) 



d'g{i, 1) 



'^^ dsedsp 



Cut 



(14) 



defines the variance between electronic and photonic 
events. The long time behaviour is given by Cg ^ q^^i 
so the correlation coefficient can then be defined as^: 



Cll 



Cr2gCr2p V(cio + 2c2o)(coi + 2co2) 



(15) 




FIG. 1: When the chemical potencial in the right lead, ^ is 
above the energy of both levels, the electron remains in the 
quantum dot and photons are spontaneously emitted analo- 
gously as photons in resonance fluorescent atoms. 



Similarly to the electronic(pliotonic) correlations, 
where the sign of the second order cunimulant, C2o(o2) i de- 
fined the sub- or super-Poissonian character of the noise, 
the sign of cn gives the character of the electron-photon 
correlations. If cn > 0, the detection of a transmitted 
electron would involve the detection of a photon in a 
short lapse of time, while cn < involves distant events. 

The electron-photon correlation coefficient is limited 
to \r\ < 1, having r = 1 for the case where the num- 
ber of detected electrons is proportional to the number 
of detected photons: Ue oc Up. r = means uncorrelated 
events. Note that independent events give r = 0, but the 
opposite is not necessarily true, as will be shown below. 
Analogously to the Fano factor for the second order cu- 
mulants, the deviation of the third cumulants from the 
Poissonian statistics can be parametrized by the coeffi- 
cient: 



Ve{p) 
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In what follows, different configurations will be dis- 
cussed concerning the relative positions of the energy lev- 
els with respect to the chemical potentials of the contacts. 
As will be shown, electron and photon fluctuations and 
their correlations are strongly sensitive to the concrete 
configuration of the system. 



III. RESONANCE FLUORESCENCE LIMIT 

The chemical potential of the left lead is considered 
to be well above the energies of the QD, so it can be 
considered infinite. If the chemical potencial of the right 
lead, /i, is also above the energies of both levels, /i > £1(2)1 
the QD is always populated by one electron and transport 
is cancelled. Then, this case is completely analogue to 
the resonance fluorescence problem, where spontaneously 
emitted photons play the role of fluorescent photons: the 
trapped electron is coherently delocalized by the driving 
fleld between the two levels performing photo-assisted 



Rabi oscillations until the emission of a photon, then the 
electron is relaxed to the lower level, cf. Fig. [TJ 

We consider a small correction to this behaviour due 
to the thermal smearing of the Fermi level for flnite tem- 
peratures. Then, there is a contribution of transport by 
a small but finite probability for the electron to be ex- 
tracted to the right lead when it occupies the upper level. 
The Fermi distribution function weighting this transi- 
tion can be approximated by /2 = a; sa g'^'^^"''', where 
/3 = {kBT)~^, see Fig. [TJ The effect of thermal smear- 
ing on electronic transport through a quantum dot has 
been controlled recentlj^. Then, photons deviate from 
the resonance fluorescence like statistics because the QD 
may be empty during short lapses of time. It would be 
the case if the resonance fluorescent atom could be even- 
tually ionized. Then, from the Taylor expansion for low 
X, one obtains a finite electronic current: 
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^'TlFr 



(72 + 2r!2) (Fl + Fr) 



Ox' 



(17) 



that introduces a small contribution in the photonic emis- 
sion: 



7f]2 



7r!2FR, (J72 + 37 (Fl + Fr)) 



72 + 2f]2 2((72 + 2r!2)2(FL+FR) 



x + O (x^) 

(18) 

The photonic resonance fluorescence behaviour as well 
as electronic transport quenching is recovered for a; = 0, 
cf. Fig. [21 There, it can be seen that the sub-Poissonian 
photon behaviour goes super-Poissonian in the vicinity 
of the resonance, as described by p^ . In those regions, 
and opposite to what is seen in resonance, the AC in- 
tensity increase the deviation of the statistics from the 
Poissonian values. From (|17p . (fT5)) and the expressions 
shown in Appendix \X\ for the second order moments, one 
obtains the contribution of the thermal smearing of the 
collector in the electronic statistics (for Fl = Fr = F) 
and the expected photonic Fano factor: 



F, 



Fr, 



1 



(Wr 



^\ (72 + 2f72) 

(72 + 2172) 



^""^ -nAx + 0{x'){19) 
2 ' Oix). (20) 



The driving field induces sub-Poissonian photonic noise 
which (in the limit x — 0) reaches a minimum Fp^m — -r 
for rim = ll ^ before the Rabi oscillations dominate 
the dynamics over relaxation processes, cf. Fig. [3] The 
electron-photon correlation coefficient becomes (see Ap- 
pendix |A| : 



r = (2F7 (72 - 101^2) - J72 (72 + 2^(2)) 



^x 



16F (74 - 2^^272 + 4f74) 



(21) 



where it is clear that the AC fleld contributes to nega- 
tive electron-photon correlations, cf. Fig. [3] The third 
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FIG. 2: /i > £2= Dependence of the photonic current, Fano 
factor and skewness with the detuning, for different field in- 
tensities in the regime where no levels are in the transport 
window: ei,2 < fJ.- Fl = Fr, = F = 1, 7 = 0.1, 2; ~ 0. Since 
electronic transport is cancelled in this regime, the photonic 
statistics are equivalent to the resonance fluorescence prob- 
lem. Sub-Poissonian photonic statistics are only found near 
resonance. It must be noted here, however, that the validity 
of these results, obtained within the rotating wave approxi- 
mation, is guarranteed only for A^^ ^ 0. 



FIG. 3: /x > £2: Dependence of the photonic current, Fano 
factor and skewness and the electron-photon correlation with 
(a) the field intensity, fl, for difi'erent photon emission rates 
and (b) the photon emission rate, 7 for difi^erent field inten- 
sities in the regime where no levels are in the transport win- 
dow: £1,2 < M- Fl = Fr = F = 1, a; « 0. Fp and rjp show a 
pronounced minimum in their dependence with the field in- 
tensity typical for resonance fluorescence. In the non driven 
case, Ip = 0, Fp = rjp = 1 and r = 0. 



cumulants become: 
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+ 0{x), (23) 



for photons. 

Interestingly, the strong photonic noise supression co- 
incides with a region where the skewness almost vanishes, 
cf. Fig. [3^, leading to the possibility to operate the de- 
vice as a regular boson source. 

Two asymptotic limits of the results presented above, 
the undriven and high field intensity limits will be con- 
sidered. 



A. Undriven case, f2 = 

In the absence of driving, once an electron occupies 
the lower level -by direct tunneling from the leads or 
by relaxation from the upper one- there is no process 
able to remove the electron from the lower level. Then, 
the stationary state of the system coincides with p2 = 
1 and both photon emission and electron tunneling are 
cancelled: 



yi,j. 



(24) 



As expected, the cancellation of photon emission makes 
all the photonic cumulants Poissonian, so Fp — rjp = I. 
However, a small contribution of the tunneling through 
the upper level modifies the electronic shot noise, so the 
Fano factor: 
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2a;rLrR 
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(Fl + Fr) (27 + xFr) - a;7FR 
and the skewness of the statistics: 

6xFlFr (2Fl (7 + xTr) + Tr (xFr - {x - 2)7)) 
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(Fl (27 + xTr) + Fr (a:FR - (x - 2)7))' 



(26) 

are spuriously super-Poissonian. Then, also the electron- 
photon correlation is affected. For the case Fl = Pr = P: 



(4 — x)^x 



2(2xP -(- (4 - x)7)(4a;P + (4 - x)-f) ' 



(27) 



All these deviations are obviously cancelled as the rate 
for extracting an electron from the upper level, xPr, goes 
to zero. 



B. High intensity limit: fi -^ 00 

Increasing the intensity of the AC field, the electron 
tends to occupy the upper level with a probability: p2 — 
Ir^Xu-xlr" ("^ 5' "^^^^n x ^ 0). Then, it can tunnel to 



the right contact (with a probabihty xFr) or be relaxed 
to the lower level (with a probability 7), contributing to 
finite electronic and photonic currents: 



/. 



2.TrLrR 



4rL + (4 - x)Tn 
7(2rL + (2 - a;)rR) 
4rL + (4 - x)Tn 



(28) 
(29) 



The electronic dynamics is then quite similar to the single 
resonant level case,i2. so the Fano factor becomes slightly 
sub-Poissonian: 



F. = 1 



SxrLFR 



(4Fl + (4 - x)TnY 



(30) 



Since the occupation probability of the upper level at 
high field intensity is maximum, so it is the probability 
of finding the QD unoccupied, po = 4rL+^4^a;)r ' after 
the extraction of the electron to the right lead. This 
introduces lapses of time when photon emission is sup- 
pressed, affecting the photonic statistics by turning it 
super-Poissonian: 



Fp = l 



2x7Fr 



(4Fl + (4 - x)Tnr 



(31) 



On contrary, electron-photon correlation is negative since 
the detection of an electron (photon) reduces the proba- 
bility of detecting a photon (electron) : when an electron 
has tunnelled out of the system (therefore the quantum 
dot is empty), photon emission is suppressed. On the 
other hand, when a photon has been emitted, the upper 
level is unoccupied and no electron can be extracted from 
the quantum dot. For Fl — Fr — F: 



/2i7(8 - 3a;) 



^(4 - a;)(64 - 24a; + x^){2x"f + (8 - x)^T) 
For the higher moments, one obtains: 
24a; (64 - (24 - a;) x) 



r]e = I- 



(8 - xY 



(32) 



(33) 



IV. 



^ , 6a;7 ((8 - a;)^F ~ (8 - 3x)7) ,„^, 

^'^ = ' + (8^xY^^ ■ ('') 



DYNAMICAL CHANNEL BLOCKADE 
REGIME 



If the chemical potential of the right lead lies between 
the energy levels of the QD, ei < /i < £2 and therefore, 
/i = l-a;, /2 = 0, wherea;«e'3(^i-^) and /3 = {kBT)-\ 
electronic transport is strongly suppressed through the 
lower level, cf. Fig. |4l Then, since only one electron is 
allowed in the system, the occupation of the lower level 
avoids the entrance of electrons from the left lead and 
the current is blocked. This mechanism, which is known 




FIG. 4: Dynamical channel blockade configuration, where the 
electronic transport is strongly suppressed through the lower 
level, though there is a small probability introduced by ther- 
mal smearing of the Fermi surface in the right lead. Again, 
the chemical potencial of the left lead is considered infinite. 
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FIG. 5: Dynamical Channel Blockade: Dependence of /e, 
Fe and r^e with (a) the field intensity, Q, for different photon 
emission rates, (b) the photon emission rate, 7 for different 
field intensities, in the dynamical channel blockade regime: 

£1 < /I < £2- M < £1.2- Fl = Fr = r = 1, n = n/r, 7 = 

7/F, X ~ 0. As discussed in the text, the super-Poissonian 
electronic Fano factor typical for dynamical channel blockade, 
is turned sub-Poissonian by the AC intensity. 



as dynamical channel blockade, predicts electronic super- 
Poissonian shot noise in multichannel systems like, for 
instance, two-level quantum dots^ or capacitively cou- 
pled double quantum dots^i^ as well as positive cross- 
correlations in three terminal devices-*' . It has been pro- 
posed as the responsible of noise enhancement measured 
experimentally in multilevel quantum dots'^^ and stacks 
of double quantum dots'^'^ . 

The blocking of the current is not forever since the 
electron in the lower level has a finite but small probabil- 
ity of tunneling to the collector, a;FR, due to the thermal 
smearing of the Fermi level. Then, the trapped electron 
eventually escapes to the right lead allowing electrons to 
tunnel through the upper level before the lower one is 
again occupied. Thus, the current is restricted to short 
lapses of time while for long periods t '^ (a;FR)^^ trans- 



fc 



0.5 
1 



-0.5, 



(a) 


',--' 


-■ ■ ■< '"^ T . 1 . 1 . 1 . — 








— ^0.1 


s . ■ 


.... ^1 . 




-- "FlO, 


- , 1 , |^~,~-1-- 


-r-T . - 




^ 




~---__^ 


- 


-.1.1.1 


^■~--___ 




- ' 1 ' 1 ' 1 


' 1 ' ■ 


\L^— "^^^ 


- 





T^r-r- 




become (see Appendix IBJ) : 







N ,.■"■' 
. ■■ ' ^ 

1 


— £2=0 

£2=1 - 
-- £2=5 




2 4 6 



10 2 4 6 

y 



8 10 



FIG. 6: Dynamical Channel Blockade: Dependence of the 
photonic current, Fano factor and skewness and the electron- 
photon correlation coefficient with (a) the field intensity, for 
different photon emission rates and (b) the photon emission 
rate, 7, for different field intensities in the dynamical channel 
blockade regime: ei < ^ < £2- A* < £1.2- Fl = Fr, = F = 1, 



port is quenched. This bunching of electrons is reflected 
in super-Poissonian shot noise. 

Photon-mediated relaxation introduces an additional 
way to occupy the lower level when current is flowing 
through the upper one, shortening the lapse of time when 
transport is allowed. Thus, the electrons are transferred 
in smaller bunches and the supcr-Poissonian character of 
the electronic noise is reduced. The detection of a pho- 
ton is always at the end of a bunch of electrons and im- 
plies the cancelation of the current, leading to a positive 
electron-photon correlation. 

The introduction of the AC field pumps the electron 
in the lower state to the upper one, giving the electron 
a finite probability to tunnel to the right lead or to be 
relaxed by the emission of one photon. This reduces 
the electronic shot noise by reducing the duration of the 
lapses of time when transport is blocked (opposite to the 
effect of photons). Thus, when a; = 0, the electronic cur- 
rent and the photonic emission are proportional to the 
driving intensity: 



2n^ri,Tj- 



FRlfl +3r!2 



Fl fR,(27 + rr 



4f^2 



Ip- 



7^(2 (2Fl + Fr) 
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rl + 3f72 



Fl FR(27 + rR; 
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(35) 



(36) 



and channel blockade is removed, see Figs. [5] and [5] We 
have defined Fr — j + Fr. Considering, for simplicity, 
the case Fl = Fr — F and F = 7 -I- F, the Fano factors 



F, 



Sn^ + 2 (272 + I5F7 -I- 9F2) rj2 - 2Ff2(37 + 2F) 



for electrons, and: 



F„ 



(tQ^ + f(37 + 2F) 



27rj2 (22F2 + 287F - f22 



F (jn^ + f (37 + 2F) 



(37) 



(38) 



for photons. As can be seen in Fig. [6^, the minimun 
that appeared in the resonance fluorescence conflgura- 
tion still appears, but its depth and position now de- 
pend on the tunneling rates. The modification of the 
resonance fluorescence behaviour is also reflected in the 
super-Poissonian large AC intensities asymptotic value 
(discussed below). 

The electron-photon correlation coefflcient will be con- 
sidered in the asymptotic non-driven and high field in- 
tensity cases. As expected, the driving contributes to 
make the electronic noise sub-Poissonian and the pho- 
tonic one super-Poissonian. However, it has to compete 
with the photon emission that contributes to bring the 
electron to the lower state and to block the current. In 
Fig. [Hb, it can be seen how the pumped electronic cur- 
rent is decreased by the photon emission rate and the 
Fano factor tends asymptotically to be Poissonian. The 
positive electron-photon correlation is decreased by the 
ac field since the emission of a photon does not imply 
transport blocking anymore, as seen in Fig. [51 



A. Undriven case, f2 = 

The most interesting features appear in the absence 
of the ac field, where the consequences of the dynami- 
cal channel blockade are maximal and there is a strong 
dependence of the statistics on the thermal smearing fac- 
tor, X. In the absence of photons (7 = 0), the electronic 
current and Fano factor are: 



le 
F, 



2xFLrR 



(x + 1)Fl + Fr 
2FL((l-a;)2FL 



(1 - 3x)Fr) 



((a; + l)FL 



rR)2 



(39) 
(40) 



It is interesting to see here how the Fano factor can be 
tuned by the asymmetric coupling to the leads: Fe — 3 (if 
Tl » Tr), Fe = 2 (ff Fl = Fr) and F, = 1 (if Fl « Fr). 
In the lastest case, the contribution of xFr is diminished 
and the left barrier controls the transport (in this limit, 
the current is le = 2a;FL). Then, the transferred elec- 
trons are uncorrelated one from the others resembling 
the behaviour of the single barrier problem briefly dis- 
cussed above. The case Fl » Fr was studied in Ref.— 
without considering the processes that introduce an elec- 
tron from the collector to the lower level, with a rate 



(1 — a;)rR. These transitions do not contribute in this 
particular hmit, but that is not the case for the rest of 
configurations. 

Considering photon emission and small x, one can ex- 
pand the first coefhcients for the electronic and photonic 
statistics, as well as for the electron-photon correlations. 
Relaxation by photons contribute to shorten the bunches 
of electrons flowing through the upper level when the 
lower one is empty, thus reducing both the electronic cur- 
rent: 
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2rLrR (7 + Tr) X 



Tr (7 + Tr) + Tl (27 + Tr) 
and Fano factor {noise reduction by noise): 

2rLrR 
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(FrF R + Tl (7 + f R 



o 



without affecting to its super-Poissonian character, cf. 
Fig. ]Ejp. Again, we defined, for simplicity, Fr = 7 + Fr. 
The photonic current and Fano factor become: 



^ _ 7rLrRX 

" ^ FR(7-hFR 



Ox' 



F„ = 1- 



Fl (27 + Frj 
27rLrR (7 + 2Fl + 2Fr) X 

(rR(7 + rR) + rL(27 + FR))' 



0^-2 



(43) 
i)44) 



while we obtain, for the electron-photon correlation co- 
efficient (in the case Fl = Fr = F): 



7(7 + F) 



(37-h2F)(37-h4F) 



0{x). 



(45) 



Interestingly, the electron-photon correlation is roughly 
independent from x, which allows to extract information 
that is not provided by the flat photonic Fano factor, cf. 
Fig. d). 

The expected positive electron-photon correlation are 
obtained. On the other hand, the presence of electronic 
transport affects the Poissonian photonic statistics by in- 
troducing a sub-Poissonian component: once a photon 
has been emitted, the electron is relaxed to the lower 
level blocking the transport. A second photon will not 
be detected until the electron tunnels to the collector and 
another one enters the upper level, so photonic events are 
well separated in time. 

From the third cumulants, one obtains: 
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6 (7rLrR (7 - 



■FL(27-f-FR))" 
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(rR(7 + rR)+FL(27 + FR))^ 
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0{x^A7) 



B. High intensity limit: £7 — » 00 

If the intensity of the driving field is large enough, the 
dynamical channel blockade is completely lifted, finding 
electronic and photonic currents: 

/. /„ 1 



2(H-a:;)FLFR 7 (2Fl + (1 - a:)FR) 4Fl + (3 - a;)FR 

(48) 
so sub-Poissonian electronic noise and super-Poissonian 
photonic noise are recovered: 
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0{x). 



(49) 
(50) 



Interestingly, in this regime, the photonic influence is 
washed out from the electronic statistics. Also, the AC 
field allows the extraction through the upper level of an 
electron that has been relaxed by the emission of one 
photon. This means that the electron-photon correlation 
becomes negative (for Fl = Fr = F): 



r = — 5i 



27 



123(27 + 49F) 



0{x). 



(51) 



Then, by tuning the driving intensity, one can ma- 
nipulate the character of the shot noise of electrons and 
photons, turning the super(sub)-Poissonian statistics to 
sub(super)-Poissonian for electrons(photons) when in- 
creasing 51. 

Higher moments are also obtained, giving: 
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V. BOTH LEVELS IN THE TRANSPORT 
WINDOW REGIME 

If the energy of both levels are above ^, £1,22 > M 
(/i = /2 = 0), the two of them contribute to electronic 
transport, cf. Fig. [7l In this particular case, quantum 
interference effects may be important^^ depending on the 
concrete geometry of the system. However, in the weak 
coupling and high frequency limit case considered here, 
£2 — El ^ Fl,r, they can be disregarded. 

Contrary to the previous regimes, the contribution of 
the empty state: 



Po = 



Tr 



2Fl + Ff 



(54) 



plays an important role here. It strongly affects the sub- 
Poissonian character of the photonic noise. 




FIG. 7: System configuration discussed in section|Vl witii tlie 
two levels in tlie transport window, £1,62 > fJ,- 



Since the tunneling rates are considered independent 
on the energy, electronic transport does not depend 
on the level that the electron occupies when tunnel- 
ing through the QD. Then, the transport characteristics 
(electronic current and noise) are independent of the field 
intensity, detunning and the spontaneous emission: 
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(2rL 



rR,)2- 



(55) 
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This case is similar to the single resonant level with a fac- 
tor 2 in the tunneling from the collector, reflecting that 
an electron in the left lead finds two different possibilities 
before tunneling into the QD. Similarly to the single res- 
onant level, the Fano factor is sub-Poissonian. However, 
the contribution of the two levels increases the noise. The 
normalized third cumulant becomes (see Appendix [Cl) : 
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i2rLrR {ATI + ri 



(2rL + T^y 



(57) 



Interestingly, the two resonant levels statistics coincides 
with the single resonant one when writing Ti^/2 for Fl. 
That is not the case for the photonic statistics, that de- 
pends on the population of the upper level and, therefore, 
on the AC field parameters. For instance, the photonic 
curret is: 



7Fl (2^2 + Fr(7 + 2Fr)) 
(2Fl + Fr) (72 + 29? + 37FR + 2F|) ■ 



(58) 



The expressions for the second order moments are quite 
lengthy, unless one considers a simpler case, where the 
tunneling rates are the same through both barriers, Fl = 
Fr = F. Then, one obtains a sub-Poissonian Fano factor: 
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which can be tuned to super-Poissonian for high enough 
intensities. The electron-photon correlation, obtained 
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FIG. 8: £i,e2 > m: Dependence of the photonic current, 
Fano factor and skewness and the electron-photon correlation 
coefficient on (a) the field intensity, f2, for different photon 
emission rates, and (b) the photon emission rate, 7 for dif- 
ferent field intensities for ^ < £1,2. Fl = Fr = F = 1. 
The electronic statistics (not shown) is sub-Poissonian and 
not affected by the AC field nor by photonic relaxation. The 
electron-photon correlation coefficient is positive if 7 > F. 



from (|C2p . may be positive or negative depending on 
the concrete parametrization of the system, as discussed 
below. In concrete, positive correlation is obtained when 
Fl ^ Fr as well as, for low intensity driving, when the 
tunneling rates are small compared to the photon emis- 
sion rate, cf. Fig. [S] 



A. Undriven case, fJ = 

The emission of a photon, in this case, depends on the 
tunneling of an electron from the left lead to the upper 
level. Then, it can tunnel to the collector directly or 
after being relaxed to the lower level by the emission of 
one photon. Therefore, photons adopt the electronic sub- 
Poissonian statistics: 



Fp = 1- 



^jVlTr (7 + 2Fl + 2Fr) 

(7 + Fr/ (2Fl + Fr)' ' 



(60) 



which is mantained for all the low AC intensities regime, 
and the resonance fluorescence-like behaviour is com- 
pletely lost, cf. Fig. M 

The sign of the electron-photon correlation depends on 
the asymmetry of the tunneling couplings. Concretely, 

in the case 7 <C Fr, it is positive if Fl < ^. Also if 
the photon emission rate is large enough compared to 
the tunneling rates. Concretely: once an electron occu- 
pies the upper level, it will rather be relaxed to the lower 
level and tunnel to the collector than directly tunnel from 



10 



the upper level. Then, the probability of detecting con- 
sequently one photon and one electron increases, thus 
making the electron-photon correlation positive if: 
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This is more clearly seen when considering Fl = Tr = F 

r == (57 - F) 



7 



10(7 + F)(772 + IO7F + 9F2) ■ 



The coefficient: 
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Vp = 
also shows sub-Poissonian behaviour 



(62) 



(63) 




B. High intensity limit: f2 ^ cx) 

For high AC field intensities, the contribution of the 
chemical potencial of the collector is only reflected in 
the occupation probabilities. Particularly important for 
the photonic dynamics is the probability of finding the 
QD in its empty and lower states, since it limits photon 
emission. The current, in this case, is: 



7rL 



2Fl + Fr 



(64) 



As seen in the previous regimes, the occupation of the 
empty state affects the sub-Poissonian statistics (ex- 
pected for resonance fluorescence) by turning it to super- 
Poissonian values: 



Fp = l 



I^R 



(2Fl + Fr) 



(65) 



Comparing to (j31[) and (j50p . the higher unoccupation of 
the QD involves a higher super-Poissonian character in 
the photonic statistics. 

High intensities allow the emission of several photons 
before the electron is extracted to the collector. Also, an 
electron tunneling from the emitter to the upper level, 
can be extracted to the collector from the lower level 
without the emission of a photon. Then, the electron- 
photon correlation tends to be negative. However, if Fl 
is small, po ~ 1 — 2Fl/Fr ^ Pi,P2, i- e., the probabil- 
ity of finding the QD empty is almost one. Then, the 
detection of photons and electrons is restricted to short 
lapses of time, which makes the electron-photon correla- 
tion positive: 
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The third order coefficient gives: 
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FIG. 9: Diagrams of the Selective tunneling configuration, 
where each level is coupled to a different lead. ( Up): the par- 
ticular system considered here. (Bottom): a possible physical 
realization by considering a triple quantum dot where the lat- 
eral ones are strongly coupled to the leads so they behave as 
zero dimensional leads. 



VI. SELECTIVE TUNNELING 
C ONFIGURATION 

A particularly interesting configuration in the high bias 
regime (/i = /2 = 0) where the electron-photon correla- 
tion is paradigmatic, needs unusual coupling to the leads: 
electrons can enter only to the upper level and tunnel 
out only from the lower one. That is: F2L = Fi^ — F, 
FiL = T2B. — 0, cf. left diagram in Fig. [9]This se- 
lective coupling to the leads could be obtained by zero- 
dimensional contacts consisting in neighbour single-level 
QDs strongly coupled to the leads^. Then, if the level 
of the left(right) dot is resonant with the upper(lower) 
level, the emitter(collector) will be uncoupled of the 
lower(upper) level, see lower diagram in Fig. [5] Any 
eventual coherence between the central dot and the lat- 
eral ones is asumed to be rapidly damped by the coupling 
to the leads. We note here that such a system can also 
be used to modulate non-Markovian dynamics by tuning 
the strength of the coupling of the lateral dots to the 
leads. 



A. Undriven case: Electron-photon identification 

In the absence of driving field, an electron that enters 
the upper level can only be transferred to the collector 
after being relaxed by the emission of one photon. There- 
fore, the electronic and photonic statistics are completely 
identical and Cio — cq^. This configuration is analogue to 
having two single level quantum dots which are incoher- 
ently coupled, giving sub-Poissonian Fano factors^ and 
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The third cumulants give, for Fl = Tr = F: 



(71) 



FIG. 10: Selective tunneling: Electronic current, Fano fac- 
tor and skewness as a function of (a) the driving intensity for 
different photon emission rates and (b) the photon emission 
rate for different field intensities. Fl = Fr = F = 1. 



Driven case 



maximal electron-photon correlation (see Appendix [D| : 

F -F - 1 27rLrR.(7 + rL + FR.) 
"^"''~ (7rR + FL(7 + rR))2 ^''^ 

r = 1. (70) 
I 



The ac field allows the tunneling of an electron to the 
collector without having previously emitted a photon as 
well as the emission of several photons from the relax- 
ation of the same electron. This makes the electronic and 
photonic currents differ, thus uncorrelating the electronic 
and photonic statistics. More interestingly, by looking at 
the dependence of the electronic and photonic currents 
with the detuning: 



FlFr (47A2 + (7 + Fr) (72 + rR7 + n^)) 



Tl (7 + Tr) (72 + 2Fr7 + 2n^ + Tl + 4A2 ) + Fr (47 A2 + (7 + Fr) (7^ + Fr7 + n^)) 
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Tl (7 + Tr) (7' + 2Fr7 -I- 2172 + F| -I- 4A2 ) + Fr (47 A^ -^ (7 + Fr) (7^ + Fr7 + r!^)) ' 
it can be seen that their second order derivative obeys: 
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(73) 
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which is reflected in a resonance to anti-resonance crossover. As a consequence, one can extract information on 
the sample-depending spontaneous photon emission rate, 7, by externally modifying the tunneling couplings to the 
collector—. The system, in this case, behaves as a photon emission rate probe. 

Considering Fl = Fr = F, for simplicity, and the resonance condition, Ac^ — 0, we obtain for the Fano factors: 
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2 (74 + 4F73 + 5F272 + 30272 + 2T^"f + 3FO27 + 20-* + 4F202 
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(75) 
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cf. Figs. [TUlandfTTl From the Fano factor, it can be seen the photons become super-Poissonian for high enough 
that the electrons obey sub-Poissonian statistics while 
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FIG. 11: Selective tunneling: Photonic current, Fano fac- 
tor and skewness as a function of (a) the driving intensity, 
for different photon emission rates, and (b) the photon emis- 
sion rate, for different field intensities. Fl = Fr = F = 1. 
Inset; Electron-photon correlation coefficient as a function 
of the field intensity for different couplings to the left con- 
tact, Fl, with Fr = 7 = 1. It is possible to tune the sign 
of the electron-photon correlation by means of the tunneling 
coupling asymmetry. 



FIG. 12: Selective tunneling: Electronic current, Fano fac- 
tor, skewness and electron-photon correlation as functions of 
the frequency detuning for different field intensities. Fl = 1, 
Fr = 0.001 and 7 = 0.001. In the insets the corresponding 
photonic values are shown. When the photon emission rate 
is much smaller than the tunneling rates and Fl > Fr, the 
system behaves as a coherently coupled double quantum dot, 
showing a sub-Poissonian minimum in the Fano factor which 
is between two super-Poissonian peaks. 



field intensities. The driving field also contributes to 
make the electron-photon correlation coefficient negative, 
see Appendix [Dl 

In the absence of relaxation, this configuration can 
be mapped into a coherently coupled single level dou- 
ble quantum dot, where interdot hopping played the role 
of the ac driving (whithin the rotating wave approxima- 
tion). Then, in the particular case where Fr < Fl, the 
noise is sub-Poissonian in resonance, having two super- 
Poissonian peaks in its vicinity, when the influence of 
photons is small, as seen in Fig. [T^l For 7 — 0^^: 
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2V?Vi^ (4 (Fr 



Fl) A2 -k Fr (2r!2 



r| 



-3r.rSP 



(rRf^2 



Fl (2^(2 + r^^ 



-4A2))^ 



This kind of features has been the subject of recent works 
in double quantum dot systems where the double peak 
structure in the electronic Fano factor becomes asym- 
metric by the effect of temperature^!^. In our case, the 
level energies are not shifted, so the contribution of pho- 
ton emission is constant all over the ac frequency tuning 
and the double peak remains symmetric. As a conse- 
quence, transport is not quenched by detuning and elec- 
tronic noise is sub-Poissonian far from resonance (where 
the ac field has no effect on transport), as expected from 
(I69p . Interestingly, the maximal electron-photon correla- 
tion observed far from resonance vanishes for A^^ = 0. 



The double peak in the electronic Fano factor is washed 
out for larger photon emission rates, even for the case 
Tl » Tr: 
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as seen in the insets of Fig. [T31 On contrary, in this 
regime, it is the photonic noise which is sub-Poissonian 
but for two super-Poissonian regions around the reso- 
nant frequency, recovering the resonance fiuorescence be- 
haviour, see Fig. [T21 



C. High intensity limit: fi — + cx3 

An intense driving involves the delocalization of the 
electron between the upper and lower level, so it has the 
same probability of occupying each of them: pi = p2 = 
rL/(2FL + Fr). In this case, the resonant currents are: 



riP 



Li R 



7rL 



2Ft 



(79) 



The electronic dynamics becomes independent from pho- 
ton emission so the Fano factor coincides with that ob- 
tained for transport through a double quantum dot in the 
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FIG. 13: Selective tunneling: Photonic current, Fano fac- 
tor, skewness and electron-photon correlation as functions of 
the frequency detuning for different field intensities. Fl ~ 
7 = 1 and Fr — 0.001. The insets show the electronic cor- 
respondents. The double peak in the electronic Fano factor 
seen in Fig. [12] is washed out by a larger photon emission 
rate. 



As discussed in previous sections, the electron-photon 
correlation is lost by the influence of the ac field. How- 
ever, if the coupling to the leads is asymmetric and 
Fr > 2rL, once the electron has tunneled out to the 
collector, the QD remains empty for a long period of 
time (compared to the lapse of time that it spends oc- 
cupied). Then, the detection of electrons and photons 
is restricted to the short periods of time, so r remains 
positive, cf. inset in Fig. [TTk : 



/^(rR,-2rL) 



V(4r£ + r^) {ATI + 4FrFl + Fr (27 + Fr)) 



The third order cumulants give: 

_ i2rLrR_(4rt + r|2 

for electrons, and: 



(82) 



(83) 



Vp = 



^ ^ 67FR (4F^ - 2 (7 - 2Fr) Fl + Fr (7 + Fr)) 

(2Fl + Fr)^ 

(84) 
for photons. 



VII. LEVEL-DEPENDENT TUNNELING 



absense of dissipation, being sub-Poissonian at resonance 
(A^ = 0): 



1 



4FLrR 

(2Fl + Tny 



(80) 



and super-Poissonian close to resonance, if Fl > Fr^iS. 
The high probability of finding the QD empty, po — 
FR/(2rL + Fr), kills the resonance fluorescence-like pho- 
ton anti-bunching and the photonic statistics become 
super-Poissonian: 



_ ^ 27rR 

''"^(2Fl+Fr)2- 



(81) 



If the left and right barriers are equal, the tunneling 
events may differ depending on which level participates. 
This can be due to the concrete orbital distribution of 
each level. Then, one has Vi^i = Vm for the couplings 
in ll]) and electronic transport can be parametrized, if 
both level are within the transport window, by the tun- 
neling rates F2 = 27rd;|V;2p and Fi = 27rd;|14ip when 
the electron tunnels to or from the upper or the lower 
level, respectively, through any barrier i^. 

The equations of motion for the generating function, 
G{t, Se, Sp) = M{se, Sp)G(t, Se, Sp), and the density ma- 
trix, p{t) = M{l,l)p{t) (after setting Se = Sp = 1), are 
then given by the matrix: 



M{Se 






r2 



i^ 
*2 





i^ 

' 2 





Ai 



eT2 \ 



* 2 


Spl 


n 


7^ 








2 


-iA^ 


'2 



(85) 



in the same matrix form chosen to write ([5]). In this case, the decoherence term is given by: A12 



r2+ri+7 



The dependence on the level which is occupied introduces the effect of the driving field and photon emission in 
the electronic current even in the high bias regime. If, for instance, Fi < F2, photon emission photon emission will 
contribute to decrease the flow of electrons. The electronic and photonic currents are, in the general case: 



/.. = 



/„ = 



(7- 



(r2 + Ti) ((r2 + Ti) (rsFi + n^) + jv^ (7 + Ti + 2r2)) 



3ri)ri 



(72 + 3f72 + 3ri (27 + Ti)) Ts + Ti (272 

7((r2 + ri)(r2ri + r!2)+7r2) 



2ri7 + 3f72) 



(7 + ari) r2 + (72 + 302 + ^^i (27 + Ti)) Ts + Ti (272 + 2ri7 + 3172) • 
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(86) 
(87) 



As expected, if Fi < r2, the emission of photons inhibits 
electronic transport. However, the opposite is not true: 
if T2 < Ti, electrons will rather tunnel through the lower 
level, thus avoiding photon emission. These two limit- 
ing cases will be further analyzed below. In the case, 
Fi ~ F2, the electronic current is independent of both 
the relaxation rate and the driving intensity, recovering 
the behaviour described in section |Vl cf. Fig. [S] 



A. Undriven case. 

In the absence of the driving field, the electronic and 
photonic currents are given by: 



/„ 



1 



Ti (7 + r2) (Fi + F2) 7riF2 27F1 + (7 + 3Fi)F2 



In this case, the difference in the tunneling rates of each 
level is enough to define the sub- or super-Poissonian 
electronic statistics: 



2F. 



F, 



(rf 



FiF2 -2F: 



(ri(7 



7r2F; 



(27Fi + (7 + 3Fi)F2)^ 



(89) 

where we have called F; = 7 -I- F^. Interestingly, in the 
absence of photon relaxation, the Fano factor depends 
linearly on the asymmetry and increases as one of the 
levels becomes uncoupled of the leads (this case will be 
considered below): 



^.„«r,.i4(£l + Il-4 



(90) 



Photon emission diminish this effect by contributing to 
make the electrons be extracted from the lower level. 
Then, the sub-Poissonian shot noise observed in the high 
bias regime, cf. Eq. ([55]) . is recovered: 



i^e(7»r.) = i 



2Fi(Fi 



(2Fi + F2 



(91) 



On the other hand, photonic statistics remain sub- 
Poissonian, independently of the configuration: 



Fr, 



1 



27F1F2 (7 + 2Fi + 2F2 
(27Fi + (7 + 3Fi)F2)' 



(92) 



The electron-photon correlation coefficient (see Appendix 
IE| shows how electrons and photons can be uncorrelated 
by the manipulation of the tunneling rates. 



B. High intensity limit: f7 ^ cx) 

As the driving field couples the two levels, it tends to 
annihilate the particular behaviour introduced by the dif- 
ferent couplings to the leads. Thus, the currents depend 
simply on their correnspondent rate: 



/, 



Fi 



1 
3' 



(93) 



while and the electronic Fano factor and skewness be- 
come independent of the tunneling couplings: i^e = f 
and 77e = ^, consistently with (|56|) . That is not the 
case for the photonic statistics, whose second and third 
moments depend on the rates: 



Fp = 1 



Vp 



27 



9(Fi 

27(7 



r2) 

9Fi 



9F2 



27 (Fi 



(94) 
(95) 



As expected, electron-photon correlation becomes nega- 
tive: 



7 



5(27 + 9(Fi+F2))' 



(96) 



C. Fi < r2 limit 

The zero-dimensional contacts introduced in the pre- 
vious section can also be employed to simulate energy- 
dependent tunneling. If both zero-dimensional contacts 
are aligned (by tuning the gate voltages of the left and 
right QDs) with the same level of the QD, transport 
through the other level will be strongly suppressed, cf. 
Fig. [141 Thus, the occupation of the off resonant level 
blocks the electronic current. 

If the levels of the surrounding QDs are aligned with 
the upper level, in the absence of driving, as soon as the 
lower level is occupied (by the relaxation of an electron 
from the upper level), transport is canceled in a high 
bias version of dynamical channel blockade. Thus, elec- 
trons flow in bunches, while photonic transport is highly 
supressed. 

Again, the driving field removes the blockade, produc- 
ing finite electronic and photonic currents: 



F2 



7 



r!2 



7r2 + 3f^2 



(97) 
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FIG. 14: Schematic diagram of the proposed setup for a level 
dependent tunneling configuration where Fi <^ F2. 



thus reducing the super-Poissonian electron noise: 

,7r2(7 + ^2? - n^j^ + 2jT2 - Tl) - 2f]4 



F, = l 



(7(7 + r2) + 3f^2)2 



which becomes sub-Poissonian for high enough driving 
intensities, cf. Fig. [121 This configuration resembles 
a single level quantum dot coupled to a localized state, 
where super-Poissonian shot noise has been predicted^. 
The obtained Fano factor recovers their result for a non- 
dissipative situation, 7 = 0. Similar models were pro- 
posed to explain enhanced shot noise in single quantum 
dots^-. For low intensities, the photonic noise is sub- 
Poissonian, resembling the resonance fluorescence but, 
for Q > ^2r2(27 + r2), the contribution of the empty 
state turns it super-Poissonian, cf. Fig. [TBI 



Fp = 1 - 2-fn^ 



2r2(27 + r2) 



n^ 



r2(72 + 3r!2 + 7r2)' 



(99) 



It is interesting to note that, though the electronic and 
photonic mean counts are proportional, their variances 
are not, which is reflected in the electron-photon corre- 
lation that gives r < 1, see (|E2|) . 

The undriven case gives a Fano factor that diverges 
when the photon emission is reduced, F^ = 1 + 2— 2-, and 



7"-i-6r^(7+r2) 
stochastic process. 



Ve 



In this case, relaxation becomes an 



D. r2 < Fi limit 

This case is similar to the previous one with the dif- 
ference that the contribution of photon emission has the 
opposite effect: the upper level is very weakly coupled 
to the leads so its population quenches the electronic 
current, cf. Fig. 1171 Therefore, relaxation by photon 
emission contributes to unblock the electronic transport. 

In the absence of driving, the electrons tend to be 
transferred through the lower level (and the system is 
reduced to the single resonant level configuration'*^), so 
there is no chance for photon emission. The introduction 
of the driving field populates the upper level thus reduc- 
ing the electronic current, cf. Fig. 1181 and giving a finite 
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FIG. 15: Level dependent tunneling: Electronic current, 
Fano factor and skewness as a function of (a) the field in- 
tensity, for different photon emission rates, and (b) the pho- 
ton emission rate for different field intensities, for the case: 
Fi = 10"^ and Fa = 1. 
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FIG. 16: Level dependent tunneling: Photonic current, 
Fano factor and skewness and electron-photon correlation co- 
efficient as a function of (a) the field intensity, for different 
photon emission rates, and (b) the photon emission rate for 
different field intensities, for the case: Fi = 10""" and F2 = 1. 



probability to photons to be emitted -thus acting as a 
photon pump, cf. Fig. 119b : 



/. = 



ri(7^ + 7ri + »^) 

272 4- 27ri + 3f72 

7r>2 

272 + 27ri + 3^2 ■ 



(100) 
(101) 



The ac field modifies the electronic Fano factor typi- 
cal from the single resonant level, F^ = 1/2, without 
changing its sub-Poissonian character but for the range 
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FIG. 17: Schematic diagram of the proposed setup for a level 
dependent tunneling configuration where r2 <^ Fi. 



Fi > V2O > 



F,,^l- 



2(7'(7 + ri)2 + f]2(3^2_r2)+2r!4) 



(272 + 2ri7 + 3n^Y 



The photonic Fano factor: 



Fp = l 



2-fn^ (72 - 5Fi7 - 



2F^ 



n^ 



Fi (272 + 2Fi7 + 3f72 



(102) 



(103) 



can be turned from sub-Poissonian to super-Poissonian 
by increasing the field intensity if Fi > i(\/33+5)7. Oth- 
erwise, it will be always super-Poissonian. The electron- 
photon correlation, calculated from (jE2[) . is always neg- 
ative. 

The third electronic cumulant varies between rye = j, 
for ri = 0, and Ve ~ 4? ^'-"' t^*^ high intensity limit, but 
it shows a deep minimum for low voltages where it is 
negative, cf. Fig. [THl The photonic one is removed by 
the AC field from T^p = 1 to the asymptotic limit: 



Vp 



27(7 -9Fi) 



27F^ 



(104) 



for rj — > 00. Then, the skewness of the photonic statistics 
can be tuned by the strength of the tunneling couplings. 
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FIG. 18: Level dependent tunneling: Electronic current, 
Fano factor and skewness as a function of (a) the field inten- 
sity, for different photon emission rates, and (b) the photon 
emission rate for different field intensities, for the case: Fi = 1 
andr2 = 10"^ 
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VIII. CONCLUSIONS 

A method for extracting the simultaneous counting 
statistics for electrons tunneling through an ac driven 
two level quantum dot and for photons emitted in the 
intradot electron relaxation processes is presented. It 
allows to calculate all the electronic and photonic cumu- 
lants as well as the correlation between fermionic and 
bosonic statistics, showing how they affect one to the 
other. For instance, photon emission is shown to re- 
duce the super-Poissonian electronic shot noise in the 
dynamical channel blockade regime. On the other hand, 
a purely quantum feature as is sub-Poissonian statistics 
in a two level photon source (resonance fluorescence) is 
lost as the electron is allowed to escape from the system. 
Our method can be applied to obtain the correlations 



FIG. 19: Level dependent tunneling: Photonic current, 
Fano factor and skewness and electron-photon correlation co- 
efficient as a function of (a) the field intensity, for different 
photon emission rates, and (b) the photon emission rate for 
different field intensities, for the case: Fi = 1 and r2 — 10~^. 



between processes of different kind affecting to the same 
system as could be spin dependent transport or three 
terminal devices. 

It is shown how the character of the electronic and 
photonic fluctuations can be manipulated by tuning the 
external parameters like the intensity of the AC field, 
the chemical potencial of the right lead or the tunnel- 
ing barriers. By this kind of measurements, information 
about electron relaxation times can be obtained. All the 
combinations of sub- and super-Poissonian noise can be 
selected in this way^. 



We present an analysis of the electron-photon corre- 
lations which gives a more complete understanding of 
the dynamical behaviour of each concrete sample con- 
figuration and the importance of relaxation processes in 
transport properties. In this sense, a configuration with 
a maximal electron-photon correlation is proposed. Ad- 
ditionally, this configuration can serve as a probe for the 
photonic emission rate. 

A triple quantum dot system is proposed in order to 
control tunneling through the central two level quan- 
tum dot, while the levels of the neightbour dots act as 
zero-dimensional leads. This way, assorted configurations 
which can be mapped to coherently or incoherently cou- 



J 
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pled double quantum dot systems, or quantum dots cou- 
pled to localized states can be achieved, providing a way 
to explore the effect of coherence in electronic transport. 
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APPENDIX A: PHOTONIC RESONANCE FLUORESCENCE 



When the chemical potential is above the energies of both levels, so the probability for an electron in the upper 
level level to tunnel to the collector is xTb., the up to second order moments are given by the coefficients (considering 
the Taylor expansion around x = 0): 



C20 



C02 



Cll 



n'Tlrl, {n' (7^ + 2n^) - 7(7 - 2n){^ + 2f}) (Tl + Tr)) 



2 ( (72 + 2r!2)3 (Tl + T^f 



3rn 



3o4 



7'^'rR 



7(72 - mn^) 



n^ 237^ - 8f72 



(72 + 2^(2)3 (72 + 2r!2)H4(7' + 2r!2)(rL + rR) 4(rL + rR)2 2(72 + 2^2^ 

-fn^TLTn (7 (72 - 10f)2) (Fl + Tr) - 02 (72 + 2^)2)) 



2(72 + 2f}2)3(rL + rR) 



-x + 0{x^) 



(Al) 

x + 0{x^) (A2) 
(A3) 



For the undrivcn and high AC intensity limits, one can give short expressions without having to do the Taylor 
expansion around x = 0. In the undrivcn case, 17 = 0, the electron-photon correlation is given by: 



Cll 



coi (xFlFr - 2cio (7 + 2rL + 2rR)) Sciocoi 



Tl (27 + xTn) + Tr {xTn - (x - 2)7) 7 + xFr 
In the opposite case, fl -^ 00, on obtains the electron-photon correlation from the coefficient: 

_ 4xcprLrR + ce ((7 - 4cp) (7 + 4rL) ~ {{x - 4)7 + Wcp) Tr) 

""" (4rL -{x- 4)Tn) (7 + xTr) 

and the skewncss of the electronic and photonic statistics from: 



C30 



C03 



64x3r3r| 



(4rL + (4 - a;)rR) 



x-t^Tr (4rL + (4 - 3x)rR) (2rL + (2 - x)Ty^ 
(4rL + (4 - x)rnf 



(A4) 



(A5) 



(A6) 
(A7) 



APPENDIX B: DYNAMICAL CHANNEL BLOCKADE REGIME. 



If the chemical potential of the collector is between the energies of each level, so the occupation of the lower 
level avoids electrons from tunneling through the upper one, until it is extracted with a rate xFr, the second order 



18 



correlations are given by (considering a; = 0): 

TlFr (7^ + 2rR7 + 40^ + r|) 



C20 = Cio- 



(Bl) 



(7 + Tr) (Fr (72 + 2rR7 + 3f^2 + rl ) + Tl (272 + 3rR7 + 4^2 + r2 )) 

2 7=^ + 2»2^ _^ i3r2 ^ + 6r^ + 8 (72 + »2) Tr + 2rL (572 + 8rR7 + 4»2 + 31 

"'''° (7 + Tr) (Fr (72 + 2rR7 + 3r!2 + r2,) + Tl (272 + 3rR7 + 4r!2 + r2,)) 

7r>2 (7 + 4rL + 3rR) - cpi (7^ + 2^27 + i3r|7 + er^ + 8 (72 + n^) Fr + 2rL (572 + 8rR7 + 4^2 + 3r2 n 
''°' " ''°' (7 + Tr) (Fr (72 + 2rR7 + 3f^2 + r2^) + Tl (272 + 3rR7 + An^ + r|)) ^*^ 

^ -7^' (7 + 4rL + 3rR) + CQiFLrR (72 + 2rR7 + 4»2 + r|) 
^" (7 + rR)(rR(72 + 2rR7 + 3r!2 + r|) + rL(272 + 3rR7 + 4r!2 + r2^)) ^ "^^ 

2coi (73 + 2^2^ + eri + (137 + GFl) r2 + 2 (572 + 4^2) Fl + 8 (72 + 2rL7 + »') Tr) 

'''° (7 + rR)(rR(72 + 2rR7 + 3r!2 + r2) + rL(272 + 3rR7 + 402 + r2)) ^ ^ 

For small x, in the undriven configuration, f7 = 0, the electron-photon correlation is given by: 

^ 7FlFr(7 + Fr)(2Fl + Fr)x ^ ^ 
(FR(7 + FR)+FL(27 + FR)f ^ ^ 

In the opposite limit, f2 — > oo, we find short expressions for the electron- photon correlation: 

2{x + 1)7FlFr(4Fl + (1 - 3a;)FR) 
- - (B6) 

and the third order moments: 

(4Fl + (3 - x)FRf ^ ' 

(x + 1)73Fr (8F2 + 2(3 - 5a;)FRFL + (3^2 - 4a; + l) F^ ) 
(4Fl + (3 - xjFr) 

APPENDIX C: BOTH LEVELS IN THE TRANSPORT WINDOW 

If the chemical potential of the collector is below the energy of both levels, the photonic shot noise and the electron- 
photon correlation can be obtained from: 

7 ((7 + 2Fr) »2 + 4Fl {n^ + Fr (7 + 2Fr))) 
^°' """^ (2Fl + Fr) (7 + 2Fr) (72 + 3Fr7 + 2^^ + 2Tl) 

2 73 + 2^2^ + 2Fr (572 + 12Fr7 + 41^2 + 8r2 ) + 2Fl (672 + m^ + 4Fr (47 + 3Fr)) 



(4Fl + (3 - a;)FR)' 




64{x + ifrlrl 




(4Fl + (3 - a;)FR)' 




(x + 1)73Fr (8F2 + 2(3 - 5a;)FRFL + (3^2 - 


-4. + l)F2} 



(CI) 



'-01" 



Cll 



(2Fl + Fr) (7 + 2Fr) (72 + 3Fr7 + 2^2 + 2F2^) 

^FlFr7 (7 + 2Fr) 



(2Fl + Fr) (7 + 2Fr) (72 + 3Fr7 + 2^(2 + 2F^) 



(C2) 



-f 2coiFlFr (57^ + 4^!^ + 4Fr (47 + 3Fr)) 

+C107 ((7 + 2Fr) n^ + 4Fl {n^ + Fr (7 + 2Fr))) 

-2C10C01 (7' + 2^*27 + 2Fr (572 + 12Fr7 + AQ^ + SVl) + 2Fl (57^ + 40^ + 4Fr (47 + 3Fr)))] 

We also obtain the skewness of the photonic statistics for the undriven case, f7 = 0: 

_ 7^FS^F| (272 + 8F2 + 7Fr (7 + Fr) + 2Fl (37 + 7Fr)) 
'" " (7 + rRy(2F. + FRy ' (^'^ 



and in the high AC intensity regime: 



7^rLFR(FR-2FL) ,„., 

Cos — K (C4) 

4(2FL-hFR)^ 



19 



APPENDIX D: SELECTIVE TUNNELING CONFIGURATION 

In the configuration describen in Sec. I VII the electronic and photonic correlations are given by: 

cio2rLrR (272 + 2rR7 + n^) 



C20 = 



C02 



(7 + Tr) (Fr (72 + rR7 + r!2) + Tl (72 + 2rR7 + 2n^ + rl)) 
_ 2 73 + 21727 + Tr (772 + 7rR7 + 4172 + r2^) + Fl (572 + 4»2 + sfr (27 + Tr)) 
""'^ (7 + Tr) (Fr (72 + rR7 + r!2) + Tl (72 + 2rR7 + 2n^ + r|)) 

C017 ((7 + Tr) r>2 + 2rL (172 + 2rR (7 + Fr))) 
(7 + Tr) (Fr (72 + rR7 + r!2) + Tl (7' + 2rR7 + 2^72 + r|)) 

2 73 + 2^727 + Tr (772 + 7rR7 + 4^*2 + r|) + Tl (572 + 4^(2 + SFr (27 + Tr)) 



"' (7 + Tr) (Fr (72 + rR7 + 172) + Tl (72 + 2rR7 + 2172 + r2^)) 

and, considering Fl = Fr = F, for simplicity: 

7 (F(7 + Tf (272 + r2) - 7(7 - 11F)F(7 + r)172 - (72 + r7 - 4F2) il^ - I76) 



Cll = 



(3f72 + (7 + r)(27 + r)r 
In the undriven case, f7 = 0, the third order coefFicents are: 



C30 



C03 = I'rtri, 



3T.3T.3 2 (7' + r^ + r^,) + 3(7Fl + 7rR + FlFr) 



(7rL + 7rR + FlFr) 



(Dl) 



(D2) 



(D3) 



(D4) 



APPENDIX E: LEVEL DEPENDENT TUNNELING 

The tunneling between the leads and the quantum dot may depend in the involved level of the quantum dot. We 
consider here the case F^l = F^r = Ti , where z = {1, 2}. In the undriven case (17 = 0), the electron-photon correlation 
is determined by the coefficient: 



Cll 



7riF2 (7 - Fl) F3 + (7 + Fi)2 r2 + Fl (272 + 3Fi7 - F?) F2 + 27 (7 - Fi) r2 



(27Fi + (7 + 3ri)r2)^ 



(El) 



We can give general expressions for the electron-photon correlation when tunneling through the one of the levels is 
supressed. For Fi <^T2. 



and, for r2 < Fi: 



Cll 



^2 7r2(7 + r2)^ - 172(72 + r2 + 67r2) - n^ 

""" " ^ (72 -f 3172 -f 7F2)3 

7172 (27FJ + (87^ + 172) r2 + (67-^ + 41727) Fl + 172 (72 + 172)) 

(272 + 2F17 + 3172)3 



(E2) 



(E3) 
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